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^ Abstract 



J> We present a minimal holographic model for s-wave superconductivity with unbalanced Fermi 

^ mixtures, in 2 + 1 dimensions at strong coupling. The breaking of a U{1)a "charge" symmetry is 

^ driven by a non-trivial profile for a charged scalar field in a charged asymptotically AdS^ black hole. 

The chemical potential imbalance is implemented by turning on the temporal component of a U{1) b 
"spin" field under which the scalar field is uncharged. Wc study the phase diagram of the model and 
comment on the eventual (non) occurrence of LOFF-like inhomogeneous superconducting phases. 
Moreover, we study "charge" and "spin" transport, implementing a holographic realization (and 
a generalization thereof to superconducting setups) of Mott's two-current model which provides 
the theoretical basis of modern spintronics. Finally we comment on possible string or M-theory 
embeddings of our model and its higher dimensional generalizations, within consistent Kaluza-Klein 
truncations and brane-anti brane setups. 
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1 Introduction and overview 



The occurrence of superconductive phases where two fermionic species contribute with un- 
balanced populations, or unbalanced chemical potentials, is an interesting possibility relevant 
both in condensed matter and in finite density QCD contexts (see [1] for a review). A chem- 
ical potential mismatch between the quarks is naturally implemented in high density QCD 
setups due to mass and charge differences between the quark species. In metallic supercon- 
ductors the imbalance can be induced by e.g. paramagnetic impurities, modeled by means 
of a Zeeman coupling of the spins of the conducting electrons with an effective external 
magnetic "exchange" field. 

At weak coupling, unbalanced Fermi mixtures are expected to develop inhomogeneous 
superconducting phases, where the Cooper pairs have non-zero total momentum. This is 
the case of the Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) phase [2]. The latter can develop, 
at small values of the temperature, provided the chemical potential mismatch 5/i is not 
too large (if the Fermi surfaces of the two species are too far apart, the system reverts to 
the normal non-superconducting phase) and not below a limiting value 6fi = 6fii found by 
Chandrasekhar and Clogston [5]. At this point, at zero temperature, the system experiences 
a first order phase transition between the standard superconducting and the LOFF phase. 
The experimental occurrence of such inhomogeneous phases is still unclear, and establishing 
their appearance in strongly-coupled unconventional systems from a theoretical point of view 
is a challenging problem. 

The imbalance of spin populations in ferromagnetic materials is also a relevant ingredient 
of spintronics, the branch of electronics concerned with the study and applications of spin 
transport. Spintronics constitutes a very important research area, due to its technological 
outcomes, e.g. in computer hard disk devices. From the theoretical point of view, the 
research in this area is vast and quickly developing. See [U |5] and references therein for 
interesting introductions to the subject. 

The roots of spintronics are based on Mott's "two-current model" [6]. Mott's model 
describes ferromagnetic materials at low temperature and in its simplest version it treats 
both charge and spin as conserved quantities, neglecting dissipative spin- flip interactions like 
spin-orbit terms. In the model, charge and spin currents flow in parallel and can be both 
induced either by turning on an electric fleld or a "spin motive force" [71 18]. The latter can 
be practically realized by a space gradient in the chemical potential imbalance, for example 
by means of a ferromagnetic-non ferromagnetic junction. The electric and spin motive forces 
can be described by means of two ("effective" in the case of the spiiij^ U{1) gauge flelds. 
This picture has been recently used in the system of degenerate free electrons subjected to 
impurity scattering in ferromagnetic conductors [10] . It is an interesting problem to see how 

^In 3 + 1 dimensions, the U{1) spin symmetry is contained in the SU{2) spin group which is a good 
approximate symmetry at low energy and temperature, see e.g. [8l |9] . 
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the two-current model can be realized in general setups (even superconducting ones) where 
a weakly coupled quasi-particle description cannot be employed. 

With the aim of providing some toy-model-based insights on these issues, we have studied 
the simplest holographic realization of strongly coupled unbalanced s-wave superconductors 
in the grand-canonical ensemble at non-zero temperature]^ Experimental evidence suggests 
that high superconductors like certain cuprates are effectively layered. Moreover it is 
expected that they develop a quantum critical point (QCP) in their phase diagram [12]. If the 
QCP displays conformal invariance, the physics at this point (and within the so called critical 
region) is effectively described by a 2+1 dimensional conformal field theory at strong coupling 
(and at finite temperature and density). A related holographic description, within a simple 
bottom-up approach, is provided by a gravitational dual model in 3 + 1 dimensions with the 
following minimal ingredients. The breaking of a U{1)a "charge" symmetry characterizing 
superconductivity is driven, on the gravity side, by a non-trivial scalar field charged under 
a U{1)a Maxwell field in an asymptotically AdS^ black hole background as in |T3l Ull [15] . 
The chemical potential mismatch, which can also be interpreted as a chemical potential for 
aU{l)B "spin" symmetry (see also p]), is accounted for in the gravity setup by turning on 
the temporal component of a Maxwell field U{1)b under which the scalar field is uncharged]^ 

The model depends on two parameters, namely the charge q of the scalar field and its 
mass m. For a particular choice of the latter (namely = —2 in units where the AdS 
radius is set to one), aimed at implementing a fermionic condensate of canonical dimension 
2j^we will show that the critical temperature below which a superconducting homogeneous 
phase develops decreases with the chemical potential mismatch, as it is expected in weakly 
coupled setups. However, the phase diagram arising from the holographic model shows 
many differences with respect to its weakly coupled counterparts. In particular there is no 
sign of a Chandrasekhar-Clogston bound at zero temperature and the superconducting-to- 
normal phase transition is always second order. This leads us to argue that a LOFF phase 
should not show up (we have checked that this is actually the case in the g ^ 1 limit). 
A different situation arises for different choices of the parameters which seem to allow for 
Chandrasekhar-Clogston-like bounds; we will not explore these choices in detail in this paper. 

Using standard holographic techniques, we also study "charge" and "spin" transport in our 
model. Essentially we turn on an external electric field as well as a "spin motive field" 
and look at the correspondingly generated "charge" and "spin" currents and J^. 
The holographic model provides a quite natural realization of Mott's two-current modelj^ 

^For an introduction to the literature on holographic studies of condensed matter systems see the reviews 
[TT] and references therein. 

•^Similarly, the holographic unbalanced p-wave setup of [H] contains SU{2) x U{1)b gauge fields. An 
C/(l)3 C SU{2) symmetry is broken by the condensation of e.g. a U{l)i C SU{2) vector J/7j which is 
uncharged under U{1)b- 

■^In 2 + 1 dimensions a weakly coupled Cooper pair has dimension 2. 

^Different holographic models with two currents have been studied in the literature. Nevertheless, as far 
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It also automatically provides a non-trivial conductivity matrix for the optical "charge", 
"spin" and other conductivities (we study heat transport, too) at zero momentum, both 
in the superconducting and in the normal phase (which can be also seen as a holographic 
"forced" ferromagnetic phase [^). The intertwining of spin and charge transport is medi- 
ated, holographically, by the interaction of the U{1)a and U{1)b Maxwell fields with the 
metric]^ This very general simple phenomenon, not related with particular assumptions 
in our model, suggests that analogous "spintronics" effects could occur in generic strongly 
coupled unbalanced models, for example in QCD at finite baryon and isospin density. 

The phenomenological bottom-up holographic approach we adopt in this paper, aims at 
providing information on some universal properties of classes of strongly coupled field theories 
with the same (broken) symmetries and scales. Details on the dual microscopic theories 
could be provided by embedding our model in full-fledged string or M-theory constructions. 
We will consider embeddings within Kaluza-Klein reductions of eleven or ten dimensional 
supergravities as well as within brane-anti brane setups. In both cases we find indications 
that a consistent top-down realization of our model requires the addition of at least another 
non-trivial real scalar in the gravity action. 



Organization of the paper and main results 

This paper is organized as follows. In Section [2] we provide a short review of the main 
features of unbalanced Fermi systems at weak coupling. In Section [3] we present our gravity 
model and write down the corresponding equations of motion with appropriate ansatze for 
the gravity fields. 

In Section [4] we describe the f/(l)^-charged Reissner-Nordstrom AdS solution correspond- 
ing to the normal non-superconducting phase, where the scalar field is zero. We discuss under 
which conditions this phase could remain stable at zero temperature under fluctuations of 
the scalar field envisaging the possibility of Chandrasekhar-Clogston-like bounds depending 



on the choices of the parameters (equation (4.14)). 



In Section [5] we present the results of the numerical analysis of the coupled differential 
equations when the scalar field has a non-trivial profile, i.e. in the superconducting phase. 
We thus describe the behavior of the condensate as a function of the temperature and the 
(Tc, 5/i) phase diagram (figure [s]), where Tc is found to be a never vanishing decreasing 
function of 6fi/fi . We also briefly comment on the (non) occurrence of LOFF phases within 

as we are aware of, the connection with the two-current model has not been explored so far. 

^In our model this is the only way in which the two Maxwell fields can couple. In the dual field theory this 
implies that the corresponding currents are only coupled through "gluon-like" loops. This feature is model 
dependent: in other holographic setups, such as those with two overlapping "flavor" D-branes, the coupling 
can happen also more directly, e.g. by means of non linear terms in the brane action [16' or by Fa A Fb 
terms [18J. The coupling via the energy-momentum tensor operator is nevertheless generically present in 
these cases too. 
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our model focusing on the large charge g ^ 1 limit. 

In Section M we study the conductivity matrix at various values of (T, Sfi) and the com- 



parison with weakly coupled spintronics. Starting from our general formulas (6.20), (6.21), 
(6.23), we provide some precise relations among the various conductivities in the normal 
phase (formula ( |6.26[ )). Interestingly, in the superconducting phase we find that in our 
strongly coupled model the DC conductivity related to the "spin" is enhanced, analogously 
to the "electric" one. Moreover, the "pseudo-gap frequency" Wgap in the superconducting 
phase is found to be non-linearly decreasing with Sfi/fi; the same happens for Ug^p/Tc. 

In Section [7] we comment on possible string or M-theory embeddings, providing evidence 
that an extra uncharged scalar is probably needed for the purpose, even in the unbalanced 
normal phase. It would be relevant to understand the physical meaning of this field and 
the role it could play within possible holographic realizations of ferromagnetic phases with 
magnon order parameters. We end up with appendices containing some review material and 
technical details. 



2 Unbalanced supercondutors 

At weak coupling, where BCS theory can be applied, the physical properties of unbalanced 
superconductors as functions of the chemical potential mismatch 5/i are well known. The 
(T, 5 fx) phase diagram is sketched in Figure [l] At zero temperature, the imbalance produces 
a separation of the Fermi surfaces of the two fermionic species, so that their condensation 
into Cooper pairs is suppressed as 6fi is increased. As it was shown by Chandrasekhar and 
Clogston p] (whose arguments we review in appendix |A]), the BCS superconducting phase 
becomes energetically unfavored whenever Sfx exceeds a limiting value Sfxi = Aq/\/2, where 
Aq is the BCS gap parameter at T = 5/1 = 0. If no other phases would develop, at Sfi = Sfii 
there would be a first order phase transition (with the gap jumping discontinuously from 
zero to Aq) from the normal to the superconducting BCS phase. 

Increasing the temperature above zero, a line of first order phase transitions emerges from 
the point (0, 6fii) in the (T, phase diagram. This line joins, at a critical point, with a line 
of second order transitions departing from (T°, 0) (at zero chemical potential mismatch, the 
phase transition is the standard BCS second order one). The critical temperature Tc below 
which the system is superconducting is a monotonically decreasing function of But this 
is not the end of the story. 

In 1964 Larkin and Ovchinnikov and independently Fulde and Ferrel [2] showed that at 
T = 0, there is a window of values Sfii < Sfi < 6fi2 for which the energetically favored 
phase is an inhomogeneous superconducting (LOFF) phase, with Cooper pairs having a 
non-zero total momentum k (see also [19] and [20] for reviews). The modulus k = \k\ is 
fixed by free energy minimization with respect to k (this amounts to setting to zero the 
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^/xi = ^ 5fX2 = 0.78Ao 



Figure 1: The (T, S/i) phase diagram for weakly coupled unbalanced superconductors. 

superconducting current) and it results to be proportional to 5//.. Its direction is chosen 
spontaneously. In general there could be crystalline phases where the fermionic condensate 

— * 

is a combination of wave functions with differently directed ki vectors. The simplest (FF and 
LO) ansatze amount to just assume a plane wave or two-plane wave (sine or cosine) form for 
the condensate. In these cases, the LOFF-to-normal phase transition at zero temperature 
is second order and Sfi2 is just slightly bigger than 6fii. In general, the width of the LOFF 
window in the phase diagram depends on the geometry of the crystal structure in momentum 
space. 

At finite temperature, a line of second order phase transitions departs from the point 
(0, (5/X2) in the phase diagram, and ends up in the previously mentioned critical point which 
thus becomes a tricritical (TCP) one. 

2.1 Unconventional superconductors 

In metallic superconductors, the Zeeman coupling with an external magnetic field (which 
induces a chemical potential imbalance) is usually negligible with respect to the orbital cou- 
pling. The latter can be naturally reduced in layered setups, such as the high Tc cuprates 
or iron pnictides, by taking the external magnetic field directed along the layers. These 
unconventional compounds are thus promising candidates for looking at LOFF-like effects, 
though their non-BCS nature makes the previously mentioned theoretical analysis not reli- 
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able. The same remark holds for unbalanced superconducting quark matter (21] whenever 
perturbation theory cannot be applied (e.g. for the estimated quark densities in the core of 
neutron stars). In this case, one can get some insight by means of e.g. Nambu-Jona-Lasinio 
(NJL)-like effective field theories, while lattice computations suffering from the so-called sign 
problem are not well suited for finite density setups. 

A promising setup where the LOFF effect can be experimentally investigated in various 
ranges of the "coupling", is that of trapped cold Fermi gases (see [22] for a general review). 
The latter experience a crossover transition between a weakly coupled BCS and a strongly 
coupled BEC (Bose- Einstein Condensate) phase as some parameter is varied. Experiments 
[23] on cold atomic gases with unbalanced populations seem to suggest that a Phase Sepa- 
ration (PS) scenario [2l| (with an homogeneous superfluid core and a normal surrounding 
shell) is realized instead of the LOFF one. However, a clear experimental evidence is still 
lacking, and the possibility of a LOFF phase emerging in some range of the coupling cannot 
be discarded. This issue has been theoretically investigated (in the mean field approxima- 
tion) in e.g. [23| by means of NJL models with four-fermion interactions. It turns out that, 
at fixed total density, when the coupling is tuned from weak to strong (i.e. in the BEC phase) 
the one-plane-wave LOFF window shrinks to a point and then disappears as the coupling is 
increased. At strong coupling, the unique homogeneous superfluid phase exhibits one gapless 
mode and has a second order transition to the normal phase. We will see that most of these 
properties are reproduced by our holographic model. 

3 Unbalanced holographic superconductors 

As we have recalled in the Introduction, in the case of metallic superconductors, a chemical 
potential mismatch between spin "up" and "down" species can be induced by the Zeeman 
effect i.e. by the interaction term "H/ = '^'-j^Hz^b'^^^ in the Hamiltonian, between e.g. 
paramagnetic impurities, modeled by an external magnetic field and the spin of the 
fermions. Here is the Bohr magneton and = diag(l,— 1). The effective chemical 
potential mismatch is given in this case by 6fi = H^^b- At low energy, the latter can be 
read as the temporal component of an effective U{1)b spin vector field i?^ (see also |9j). The 
two fermionic species have opposite "charges" with respect to U{1)b, hence the up-down 
Cooper condensate (charged under an electromagnetic U{1)a vector field whose temporal 
component effectively provides the mean chemical potential /x) has zero total spin charge. 
This simple picture lies at the basis of our effective holographic model and can be applied 
to more general unbalanced setups (e.g. to QCD-like ones). 

The simplest gravity model in 3 + 1 dimensions aiming at implementing holographically 
the main features of the quantum critical region of s-wave unbalanced unconventional super- 
conductors in 2 -|- 1 dimensions is thus described by the following action (see appendix[B|for 



7 



a generalization of this setting to any d > 2) 



S = 7^ / dx^y/^ 



(3.i: 



which is a very simple extension of that introduced in [HI [15] for the balanced case. The 
Maxwell field Aa (resp. Ba) with field strength F = dA (resp. Y = dB) is the holographic 
dual of the U{1)a "charge" (resp. U{1)b "spin") current of the 2 + 1 dimensional field theory; 
when y = the system reduces to that introduced in [HI [15] . The metric is mapped into 
the field theory stress tensor. Finally, the complex scalar field if) is dual to the condensate. 

The action is chosen so that it admits an AdS^ solution of radius L when all the matter 
fields are turned off. Since we are interested in finite temperature setups we will focus on 
asymptotically AdS^^ black hole solutions. Notice that the fields in the action are dimension- 
less and the gauge couplings have been reabsorbed in the overall gravitational constant k\. 
Hence the U{1)a charge q of the scalar field has dimension of an energy. 

The functional form of the potential ^(I'j/'l) is not apriori constrained by the symmetries 
in the model. For simplicity, as in [HI [15] , we will consider a potential containing only the 
mass term 

V{\il}\) = m^^^. (3.2) 

By means of standard AdS/C¥1 map, the scalar field %() results to be dual to a charged 
scalar operator of dimension 

A(A-3)=m2L^ (3.3) 

With the aim of describing a fermionic Cooper pair-like condensate, Oa = ^^'^ with canon- 
ical dimension A = 2 , we will mainly focus on a particular choice for the mass of the scalar 
field: m? = — 2/L^. Of course there is no reason to believe that this will be actually the 
dimension of the order parameter driving superconductivity at strong coupling, so we will 
also consider other mass values at T = 0. In any case the mass parameter will be taken to 
be above the the Breitenlohner-Freedman (BF) bound in AdS^, 

m^L^>-^. (3.4) 
- 4 

3.1 Ansatz and equations of motion 



From the action (3.1) one obtains the following equations of motion: Einstein's equations 



Tj 9abR "igab c-\ 
riab ^ JY--^J^ab, lO-Oj 



8 



where 



^ab = -FacF^-YacY.'+^-gabF.dF'^'+^-gabYcdY"' 



+gaiy{\ij\) + gah\dij - iqAij\^ 

-[{da^l^ - iqAaij){d,iP^ + iqA,iP^) + (a O b)] , (3.6) 



Maxwell's equations for the Aa field 



V-9 

the scalar equation 



^ :a„(v^rt"^6c) = igg'^iip^dci^ - iQ^c^) - i^idci^^ + «?^cV'^)] , (3.7) 



a„[x/=^(afe^ - iqA,xl^)g^'] + ig5"M,(a„^ - iqA^^P) + ;^rn V^'(|^|) = , (3. 



and Maxwell's equations for the field 

1 



5a(^r/^nc) = 0. (3.9) 



Let us now look for static asymptotically AdS black hole solutions of the previous equations. 
These solutions will be dual to the the equilibrium phases of the dual quantum field theory. 
For our purposes the most general ansatz for the spacetime metric is 

ds" = -g{r)e-^^'-^de + —{dx'' + dy^) + ^ , (3.10) 

together with a homogeneous ansatz for the fields 

i) = ipir), Aadx"" = (j){r)dt, Badx" = v{r)dt . (3.11) 

We will focus on black hole solutions, with a horizon at r — rn where g{rH) — 0. The 
temperature of such black holes is given by 

^^g_yrH)^ _ .3_^2) 

Using one of Maxwell's equations we can safely choose '4> to be real. The scalar equation 

bGCOHlGS 

g r 2] 2g g"^ 
Maxwell's equations for the field become 



^" + ^,'('L + i_:^)_^ + ^l-q^=0, (3.13) 



the independent component of Einstein's equations yield 

1 e^(£^+^ 1 3 y(^) eW0! . ..... 

X' + r^'^ + r ^ ; ^ = , (3.16) 
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Maxwell's equation for the v field becomes 

^" + t;'0 + |)=O. (3.17) 

When v{r) = these equations reduce to those found in [IS]. As already mentioned we will 
specialize to the case where the scalar potential contains only the mass term. In appendix [B] 
we have reported the same equations of motion for a general dimension d of the spacetime. 
In the following we will work in units L = 1, 2k,1 = 1. 



3.2 Boundary conditions 



In order to find the solution to equations (3.13 3.17) one must impose two suitable boundary 



conditions: one in the interior of the spacetime at r = and one at the conformal boundary 
r = oo, where we require AdS asymptotics. The analysis here is the standard one [T11[T2]. At 
the horizon, both g{r) and the temporal components of the gauge fields should be vanishing 
|14j . Hence we will require 

4>{^h) = ^{rH) = girff) = , and '>p{j'H),xij'H) constants. (3.18) 

The series expansions of the fields out of the horizon th, implementing the above boundary 
conditions, are the following 

(pH{r) = 4>Hi{r - rn) + 4>H2{r - rnf + • • . , (3.19) 

ipH{r) = i)Ho + i^Hi{r - th) + ipH2{r - rnf + . . . , (3.20) 

Xuir) = xm + Xmir - th) + XH2{r - ruf + • • • , (3.21) 

9H{r) = gmir - th) + gH2{'r - rnf + (3.22) 

vnir) = VHi{r - Th) + VH2{r - rnf + .... (3.23) 

At the conformal boundary we must impose a leading behavior according to the correspond- 
ing dual boundary operators. Choosing m?L'^ = —2, the scalar field should approach the 
boundary in the following way 

C C 

^jJ(r) = — + ^ + ..., as r ^ oo. (3.24) 
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With a homogeneous anstaz (3.11) Ci and C2 are constants, independent on the field theory 
coordinates x^. Our choice of mass does not lead to non-normalizable modes (— | < m?L'^ < 
— |), hence we can in principle choose whether the leading or the sub leading behavior in 



(3.24) should be the source of the dual operator O. Since we want the condensate to arise 



spontaneously, we shall require either one or the other independent parameter in (3.24) to 
vanish. Specifically we will choose 

Ci=0, (0) = v^C2, (3.25) 

where the factor \/2 is a convenient normalization as in [13] . 
Vector fields at the boundary are given by 

0(r) = /i — - + ... as r— i-oo, (3.26) 
r 

S p 

vir) = 6fi h... as r— t-oo, (3.27) 

r 

where /i (resp. p) and Sfi (resp. 6p) are the mean chemical potential (resp. the mean charge 
density) and the chemical potential mismatch (resp. the charge density mismatch) of the 
dual field theory]^ The reason behind the latter identifications has been explained before: 
a chemical potential mismatch can be realized by turning on a chemical potential for an 
effective U{l)s "spin" field, under which a Cooper-like order parameter (whose gravity dual 
is the scalar field ip) is uncharged. The U{1)b Maxwell field (of which v{r) is the electric 
component) is the holographic realization of such a field. Of course our gravity model just 
provides an effective description of the symmetries and the order parameters of the dual field 
theory. As such, the U{1)a and U{1)b fields that we treat as the holographic duals of the 
"charge" and the "spin" currents respectively, could be actually mapped into any couple of 
abelian global symmetries of which one can be broken by a vev of a charged scalar, while 
the other stays unbroken. 

Notice that we work in units L = 2k\ = 1, where the bulk fields Aa, Ba and the parameters 
/i, have mass dimension 1, while if) is dimensionless; p and 5p are charges per unit volume 
in the (2+l)-dimensional field theory, hence have dimension the radial coordinate r has 
dimension 1 in mass. 

The functions in the metric ansatz should have AdS^ asymptotics as in [T5] 

gir) = \- . . . as r ^ 00 (3.28) 

2r 

X{r) = + ... as r ^00, (3.29) 
where e is holographically mapped to the energy density of the dual field theory. 



■^In this paper we shall work in the grand-canonical ensemble with fixed chemical potentials. 
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4 The Normal Phase 



A simple solution to the equations of motion (3.13 3.17) corresponds to the normal phase in 
the dual field theory. This is characterized by a vanishing vacuum expectation value of the 
condensate O, corresponding to a vanishing scalar field '0 = in the bulk. The corresponding 
background is that of a f7(l) ^-charged Reissner-Nordstrom-^(i5'4 black hole, with metric 

ds'' = -f{r)dt' + r\dx' + dy^) + j^, (4.1) 

Here th is the coordinate of the black hole outer horizon. The gauge fields are given by 

0(r) = Jl-'^)=^,-P, (4.3) 
\ r J r 

v{r) = 5/ifl-^^ =5/i-^. (4.4) 

The temperature reads 
from which we get 



'"^12-^^^), (4.5) 



3 6 

The Gibbs free energy density reads 



th = ^TiT + I v^l67r2T2 + 3(/i2 + 6fi^) . (4.6) 



2^ 



^n = -rj,(l+ ^" (4.7) 



Notice that, due to formula (4.6), this is a function of T, /i and 5jj. 

The doubly charged AdS black hole solution is certainly not new in the holographic con- 
densed matter literature. For example, it also describes the normal phase of unbalanced 
p-wave superconductors As it was noticed in P, the normal phase of a model like 

ours can be also seen as a rough holographic realization of a "forced" ferromagnet, where 
the "spin" density Sp is supported by a non zero value of 6fi (and so, equivalently, by an 
external magnetic field). Indeed in our case (5p = if = 0. g 

At T = the doubly charged RN- 74^5*4 black hole becomes extremal and (see equation 
(4.5)) the horizon radius reads 

rl = ^{Sfi' + ft') at T = 0. (4.8) 



^In a real ferromagnet, instead, the spin density Sp is spontaneously generated. Moreover there is a non 
vanishing ferromagnetic order parameter (a magnon). The holographic description of such a setup is an 
interesting open issue. 
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In the near-horizon hmit the metric reduces to that of an AdS2 x i?^ background with AdS2 
radius given by L^2) = L'^/Q (see also Appendix b1) 
The charge density imbalance at T = reads 



ip=^q^i,. (4.9) 

Notice that this is zero at 5/i = as it happens for the normal phase at weak coupling (see 
appendix [a]) . The susceptibility imbalance ("magnetic" susceptibility) reads thus 



Therefore, in the limit Sfi <^ /i, the Gibbs free energy density of the normal phase at T = 
goes at leading order as 

a;(<5/i)^a;(0)-^-^5/i^ (4.11) 



Following the same reasonings as in [3] (see also appendix |A]) , we can argue that, provided 
a superconducting phase exists at T = 0, and it has 6p = for every 6fi, a Chandrasekhar- 
Clogston bound at T = should naturally arise also within the holographic setup. 

Let us now ask whether there are conditions under which, lowering the temperature, a 
superconducting phase {ip 0) might arise with a formation of a charged condensate below 
a certain critical temperature Tc. 



4.1 A criterion for instability 



In our model we can find a simple condition (see also [9l [26]) on the external parameters in 
order for the normal phase to become unstable at T = 0. Let us consider a fluctuation of the 
complex scalar field ip, charged under U{1)a, around the extremal [/(I) ^-charged RN-AdS* 
background. Its equation of motion has the form given in ( 3.13[ ) with background metric 
given in (4.1) and (j){r) given in (4.3). The horizon radius is fixed as in (4.8). 

In the near-horizon limit it is easy (see j26 



to see that the equation for ip reduces to that 
2g2 



of a scalar field of mass 



m 



efr(2) 



m 



1 + ^-4 



(4.12) 



on an 74^5*2 background of squared radius 



"(2) — 1/6- The instability of the normal phase 
in this limit, is thus mapped into the requirement that the above effective mass is below the 
AdS2 BF bound 

^ ^ <-^, (4.13) 



7-2 2 
-^(2)"^eff{2) 



6 



"^eff(2) 
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f){'n' + l]<2f. (4.14) 



which leads to 



For an analogous formula in the general d + 1-dimensional case see appendix |B] (see also [9J). 

When (m^ + f) < O5 i-^- ""^^ < — |, the instability occurs for every value of This 
will indeed be the case for = —2. This suggests, quite surprisingly, that in these cases a 
superconducting phase with non-trivial scalar profile could be always preferred at T = 0, no 
matter how large is the chemical potential mismatch]^ 

In the cases in which (m^ + |) > 0, instead, the normal phase will show instability when 

which gives an actual bound on 6fi/^ provided that > 2m^ + 3. The above condition 
resembles the Chandrasekhar-Clogston bound of weakly interacting superconductors. 

According to the comments in [27] , we expect that a violation of the 74^5*2 BF bound leads 
to a continuous phase transition of the Berezinskii-Kosterlits-Thouless (BKT) type at T = 
(see also [28]). In BKT transitions the order parameter goes to zero exponentially instead 
as with the power law behavior of second order phase transitions. Around these phase tran- 
sitions there should be a turning point in the phase diagram, with the critical temperature 
slowly going to zero as an external parameter (for us Sfi/fi) is increased. Actually, a BKT 
transition should become of second order at T > 0. Moreover, in [28] it has been observed 
that, in a holographic model, a BKT transition at T = can only occur when the theory 
has two control parameters with the same dimension. This is precisely what happens in our 
case. Finally, notice that if the normal-to-superconducting phase transition at T = is a 
continuous one (e.g. a BKT one) the critical value of the parameter as deduced from 
the BF bound in AdS2 should correspond to the critical value at which the phase transition 
occurs (see also analogous comments in [16]). 



5 The Superconducting Phase 

If the normal phase becomes unstable at low T, we must search for another static solution 



to the equations of motion (3.13 3.17) where the scalar field is non-zero. In the dual field 



theory this corresponds to turning on a vacuum expectation value of the condensate leading 
to a spontaneous symmetry breaking of an electromagnetic symmetry and the consequent 
emergence of a superconducting phase. 



similar result was found in [9 studying the instability of an extremal dyonic black hole, electrically 
and magnetically charged under a U{1) field. 
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In the following we will discuss the results of a standard numerical analysis of the full set 
of equations of motion subjected to the boundary conditions discussed in Section 3^ The 
analysis is mainly based on the shooting technique and it is strictly valid at T > 0. We defer 
the study of the T = case, along the lines considered in [26], to future work. 



5.1 The condensate 

Let us concentrate on what we can learn from the numeric solutions. First of all let us find 
an expression for the temperature as a function of the horizon values of the various fields. 



From the general expression (3.12) and the Einstein equation at the horizon, it follows that, 
for generic 



167r 



:i2-2mVio)e ^ - -^e^(0^i + 



H 



(5.1) 



The critical temperature is found by setting (O) ~ C2 = 0|^ 

Our numerical analysis gives rise to the following results. For small values of the chemical 
potential mismatch 5^ = 0.01 and for different values of the external parameter q we obtain 
results similar to [IS]. A condensate arises below a certain critical temperature signaling 
a phase transition from a normal to a superconducting phase. The general form of these 
curves is similar to the ones in BCS theory, typical of mean field theories and second order 
phase transitions. The value of the condensate depends on the charge of the bulk field q. 
However, as in [15] , it is difficult to get the numerics reliably down to very low temperatures. 

Allowing for higher values of the chemical potential mismatch Sfi, we obtain analogous 
behaviors for the condensate. Increasing the value of 5/i (see Figure [2]) we obtain a decreas- 
ing value of the critical temperature. The phase transition is always second order. The 
most interesting result is the plot of the critical temperature normalized to T° (the critical 
temperature at zero chemical potential mismatch), against 6fi/fi. The second order phase 
transition at zero chemical potential mismatch develops inside the (Tc, Sfi) phase diagram. 
As it is shown in Figure [sj the critical temperature decreases with a qualitative fea- 

ture which we have seen also in the weakly coupled case. However, differently from the 
weakly coupled case, there is no finite value of for which Tc = 0|^ Hence, there is 

no sign of a Chandrasekhar-Clogston bound at zero temperature. This result matches with 



the expectations coming from the formula (4.14), which actually suggested the absence of 



such a bound for = —2. However, it should be desirable to refine our numerics around 
T = as done in [2Bj to definitely confirm this conclusion. In any case, we believe that it is 
unlikely that the curve in Figure [3] will suddenly drop to zero with another fiex. The phase 



dimensionless combination involving T and the chemical potentials is T/{jj^ + (5/i^)^/^. 
^^A similar phase diagram appears in [29^ in holographic superconductors in the presence of an external 



magnetic field. 
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Figure 2: The value of the condensate as a function of the temperature at /i = 1, g = 2. 
From right to left we have 5// = 0, 1, 1.5 and the critical temperature is decreasing. 



transition we find is always second order as we have also checked by a standard holographic 
computation of the free energy along the same lines as in [15] |^ Together with the absence 
of a Chandrasekhar-Clogston bound, this leads us to argue that a LOFF phase is unlikely 
to develop. 



5.2 Comments on the LOFF phase 

Let us see whether the occurrence of inhomogeneous phases is actually forbidden in our 
model. For simplicity, we will limit our analysis to the so-called "probe" approximation as 
first considered in It consists in rescaling the scalar as well as the Aa gauge field by 
the charge 1/q and taking the g S> 1 limit. In this way the backreaction of these fields on 
the metric can be neglected. This limit could only capture the physics at temperatures close 
to the critical temperature Tc, where the condensate (hence the field tp) is actually small. 
In particular, it could be reliable around a tricritical point, in case this is displayed by the 
(T, 6fi) phase diagram. 

Looking for LOFF-like phases in this limit corresponds to looking for gravity solutions (on 
a fixed background) where the scalar field spontaneously acquires a dependence on e.g. one 
of the spatial directions. In particular we will focus on complex one-plane wave solutions 



and on real sinusoidal ones, In our case the fixed background is given by a [/(l)^- charged 



^^Notice that, as we have previously argued, the phase diagram could change for different choices of the 
mass parameters in our model. We defer this analysis to future works. 
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Inhomogeneous phases of this kind in holographic p-wave superconductors have been considered in, e.g. 
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5ni ^x 

Figure 3: Second order phase transition line in the (T^, plane with q = 1, fi = 1.87. 
There are always values of below which a superconducting phase arises. 
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Let us now consider a simple single plane wave inhomogeneous ansatz 



-ixk 



A = At{r)dt + A^{r)dt . 



(5.2) 

(5.3) 
(5.4) 

(5.5) 



Inserting the above ansatz in the equations of motion (3.7) and (3.8) on the fixed U{1)b- 
charged RN-AfiS" background (5.2), (5.3), (5.4), one gets Maxwell's equations 



diAt + -drAt 

r f 

2^2 



- A = , 



and the scalar equation 



/ 



{k + A,)=0. 



+ ^ - 



f f 



ik + A.J^ 



0. 



(5.6) 
(5.7) 

(5.8) 



EOlfTHl 
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These equations always admit a trivial solution with = —k which is just related to the 
k = one by bulk gauge transformations. Both solutions thus correspond to the same 
zero-current homogeneous phase. In order to find non-trivial solutions corresponding to an 



absolute equilibrium state with zero superfiuid current we could require the Maxwell field 
Ax{r) to have UV boundary conditions such that 

AJr oo) ^ -k+ - , with J = 0. (5.9) 
r 

The large-r asymptotics for and At are taken as in the homogeneous case. From the 
above equations it is easy to realize that non-trivial solutions for A^ satisfying the boundary 



condition (5.9) are not admitted. This is true for generic values of m^, for which the scalar 
field \E'(r), dual to the same kind of operator we have considered above, has to go like 
\1' ^ Cr~'^ at large r (here A = (3 + \/9 + 4m2)/2 is the operator UV dimension). The 
same conclusion is reached starting from a more general ansatz in which At = At{x,r) and 



Ax = Ax{x, r), with the same large r asymptotics as in (5.9 ). Maxwell's equations in this case 
imply separability dxdrA^ = and consistency forces A^ and At to loose their dependence 
on X, thus reducing the setup to the one above. 

It is also possible to show that two- wave real condensates going like ip = "^{r) cos{kx) are 
excluded already at the level of the equations of motion. 

All in all, this analysis seems to exclude the possibility that a LOFF-like phase, and thus 
a related tricritical point in the (T, 6fi) phase diagram, can be displayed by our model. 



6 Conductivities: holographic spintronics 

In this section we present the study of the conductivities of the system as functions of the 
frequency of the applied external field perturbations. In this paper we limit the analysis to 
zero-momentum perturbations. Thanks to the rotational invariance of the theory in the x — y 
plane, it is sufficient to consider the conductivities in one direction, say along x. According to 
the AdS /CFT dictionary, the calculation is performed through the study of the response of 
the gravity system to vector perturbations in the x direction. The holographic computation. 



which is fairly standard, is described in Section 6.2 the un-interested reader can safely skip 



it (apart from formulas (6.20), (6.21), (6.23)) and go to the results in Sections 6.3 and 6.4 



6.1 The conductivity matrix 

The system at hands includes two U{1) vector fields describing two currents of the dual 
boundary theory. As mentioned in the Introduction, it represents a minimal holographic 



^^This is a necessary condition for having a LOFF ground state [2] and it marks a difference w.r.t. the 
standard superfiuid phases with non zero superfiuid velocity, see e.g. [STj I32j. 
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description of the "two-current model" [SI El E] which hes at the roots of the theoretical 
study of spintronics. The two-current model is based on the observation by Mott |6] that at 
low temperature the two currents of spin "up" and "down" electrons can be considered as 
two separate entities, with two corresponding conductivities which are in principle different. 
The system is typically a ferromagnetic metal with a non zero spin polarization 6p. 

In such case, an applied external electric field E'^ (providing an "electric motive force") 
does not generate only an electric current with conductivity ace-, but can also generate a 
net spin current J**. The corresponding conductivity is sometimes called (somehow asym- 
metrically) the "spin conductivity" cTcs- Moreover, this effect is reciprocal: a "spin motive 
force" generated by an external applied field E" (essentially, a dynamical gradient of popu- 



lation imbalance, V(5/i), on top of creating a spin current with conductivity ass (we call 
this the "spin-spin conductivity"), induces an electric current J'^ as well, with conductivity 
Csc [H]; the latter is precisely equal to the spin conductivity acs in time-reversal invariant 
settings. We will not consider dissipative effects bringing to spin-relaxation in this paper. 

The electric and spin motive forces can be described by means of two U{1) gauge fields: 
the electro-magnetic one and the spin one; the latter is an effective gauge boson. In this 
paper we have presented a gravity theory which precisely encodes the minimal ingredients of 
a macroscopic description of the two-current model: a charged environment, i.e. the charged 
black hole, with two vector operators, dual to the two gauge fields Aa and Ba- The two gauge 
fields are not directly coupled in the gravity Lagrangian, but their fluctuations are coupled 
via their coupling with the metric. As a result, the two dual currents are coupled through 
their mixing with the energy-momentum tensor operator of the field theory. Obviously 
enough, there is a mixing with the heat current as well. 

Thus, we would like to stress the fact that from the dual gravitational point of view, 
the mixing of the two current operators, which causes the crucial existence of the "spin 
conductivity" acs = o'sc, is completely straightforward and universal. In fact, both currents 
bring some momentum, thus they source the momentum operator Tt^; the latter is dual 
to the metric component gtx which, being a vector perturbation of the metric, mixes with 
both gravity vector fields, coupling the dual operators. In other words, this phenomenon is 
not very sensitive to the details of the gravity Lagrangian used to describe the two-current 
model: generally, a gravity theory with two conserved f/(l)'s will provide a non- zero dual 
"spin conductivity" . 

Furthermore, the gravity theory at hands includes a charged operator (under the "electric 
charge") which condenses at sufficiently small temperature. Thus, the system describes both 



^^This can be generated by "dynamical magnetic textures" , e.g. dynamical magnetizations in ferromag- 
netic conductors. Typical potential differences are generated by appropriately engineered sequences of layers 
of materials with different magnetic properties. Thus, the "spin motive force" is generated perpendicularly 
to the layers. In our gravity dual model, we instead consider the "spin motive force" in the plane; for our 
purposes the difference of the two settings is irrelevant. 
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the normal phase of a two-current model, and a superconducting phase thereof. It must also 
be remembered that the gravity theory provides a dual description of strongly interacting 
microscopic degrees of freedom, as opposed to the usual weakly interacting fermion systems 
of the standard spintronics literature. 

Actually, the gravity description does not rely on the specific microscopic origin of the two 
f/(l) currents in terms of charge and spin, although it describes the same basic features. In 
this sense, it is more "universal" and could be possibly applied to other microscopic theories, 
such as QCD-like ones. Thus, we chose to call , E"^ the current, external field and 
conductivity associated to the U{1)a under which the scalar operator is charged; these are 
the quantities which we called "electric" above (i.e. electric current J'^, field E"^, conductivity 
(Tec). Analogously, we call J^,E^,aB the quantities associated to the U{1)b under which 
the scalar operator is un-charged; above we referred to them in relation to the "spin" (i.e. 
they are J^, E^, ass in the case of electron spin unbalance). Finally, we name 7 the mixed 
conductivity (fXcs = Csc above). Nevertheless, in the discussion below we will often indulge 
with the "electric/spin" terminology for a (hopefully) clearer presentation. 

According to the discussion above, all the conductivities are included in the general non- 
diagonal matrix (see e.g. [S] and [U] for the notation) 

' (7 a aT 7 

aT kT /3T I • I I . (6.1) 




7 f^T aB 

In this matrix, (7^,0"^ are the diagonal conductivities associated to the two t/(l)'s. In 
the language used above, the former is the standard electric conductivity, measuring the 
proportionality between an applied external electric field E^ and a generated electric current 
J"^. Analogously, the "spin-spin" conductivity is the proportionality coefficient between 
a gradient in the imbalance chemical potential 6fi (which we called E^ above to underline 
the similarity with E^), and the generated spin current . The third diagonal entry in 



(6.1), K, is the thermal conductivity, i.e. the proportionality between the thermal gradient 



— ^ and the heat current Q = Ttx — fiJ^ — SfiJ^. 

Moreover, a and f3 are the thermo-electric and "thermo-spin" conductivities:^ they are 
associated to the transport of heat in the presence of an electric, or spin, potential even 
without a temperature gradient VT. Finally, 7 measures e.g. the spin current generated 
by an applied electric potential E^ even without an applied external field E^. 



The matrix (6.1) is symmetric because of time-reversal symmetry [8| ITS]. 



is sometimes called the "thermo-magnetic" conductivity |8j. 
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6.2 Holographic calculation of the conductivities 



We follow closely the presentation in [TU [15] throughout all the present section, underlining 
the novelties of the unbalanced case. In the holographic approach, the fluctuations of the 
fields Ax, Bx, gtx at the AdS boundary, act as sources for the currents J^, and the stress 
energy tensor component Ttx- Conductivity is a transport phenomenon, hence it requires 
a real time description. Since the black hole solutions we consider are classical, we must 
require in-going boundary conditions for the fields Ax, B^, gtx at the horizon. Let us take a 
simple e"*"^* time dependence for the fluctuations and consider the related linearized Einstein 
and Maxwell equations on the background 



A''.+ 



2 



A' + 



2g2^2 



A^ 



-g'tx + -9tx 
r 



B:+{1-^\B', + ^e^Bx = -e^ I -g^x + -gtx 
g 2 / 0^ Q \ r 



g'tx - -9tx + (t^'Ax + v'Bx = . 
r 



(6.2) 



(6.3) 



(6.4) 



The prime represents the derivative with respect to the bulk radial coordinate r. Substituting 

2gV~ 



(6.4) into (6.2) and (6.3) we obtain 



9 



9' 



9 



Ax - -e^ {BxV + Axcj)') = , 
9 



(6.5) 



B'' + 



2 



Bl + ^e^Bx - -e^ {Bxv' + A^ 



0. 



(6.6) 



9_ 

9 J ^ 9^ 9 
In this way we can deal with two equations in which the metric fluctuations do not appear. 
Notice that the substitution has lead to a system of linear differential equations in which the 
two gauge fields Ax and Bx are mixed. It is important to underline the role of the metric in 
such mixing: indeed, in the probe approximation, no coupling between the different gauge 
fields occurs. 



In order to solve (6.5) and (6.6), we assume the following near-horizon behavior ansatz for 



the fluctuation functioni 



17 



Axir) 



Bxir) 



r 

\ iauj 

r . 



ao + ai ( 1 ^ j + 



bo + hil 



r 
rH 



(6.7) 
(6.8) 



^^It can be checked that in order to have a non-trivial solution of the equations around the horizon, the 
frequencies of the two modes must be equal. Moreover, inspection of the behavior of the equations near the 
horizon dictates that the exponential coefficients "a" of the leading terms of the two modes must be equal 
too. 
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The IR solution depends on two integration constants, oq, 6o- The overall scaling symmetry 
{Ax — )■ XAx,Bx — 7- XBx) of equations (6.5), (6.6) allows one to fix e.g. ao = 1. 
In the UV, the behavior of the fields is 





= A^ 4 


-1^ + .:. 


(6.9) 




= Bi'^ 4 




(6.10) 


9tx{r) 




~ rdtx + ••• • 


(6.11) 



With this notation, the solution for gtx in equation (6.4) can be expressed as 



gtx 



b'Ax + v'Bx 



so that 



9 to: 



Moreover, from linearity and symmetries of (|6.5|) and (|6.6|) it follows that, on-shell, 



5?) 



iio-^Af + iuj(TbB^^^ 



(6.12) 



(6.13) 



(6.14) 



where the reason behind the labeling of the cj-dependent coefficients will be clear in a mo- 
ment. 

The linear response of a current 3°^ to perturbations (p^J^ driven by external sources 0^ 
is given in terms of retarded correlators. Formally, in momentum space, (J") = G'j?[J'^, J^](t)h- 
In our case, the retarded correlators Gr, which are proportional to the conductivities, can be 
holographically computed using the on-shell gravity action for the dual fields as generating 
functional. The on-shell action for the quadratic fiuctuations of A^^B^^ gtx can be expressed 
as just a boundary term at infinity and, after performing the holographic renormalization 
procedure, it can be reduced to 



Squad — j d^x (^A^^^ A'^P + -Bf^B^^^ — Sg^J g^J — -gtx^gtJ^ , (6.15) 

with A'^\Bi^\gil^ given in ( |6JL3| and ^I^ . 

The conductivity matrix can be thus deduced using the holographic relations 

^ Squad 



Q = 



^ Squad 



^ Squad 



(6.16) 
(6.17) 
(6.18) 
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provided the following formulas (see Hartnoll's and Herzog's reviews in [TJ for details) are 
employed 

Ef = M4°^+rf), i^f = M5i°^+rf), -^ = ^rf- (6.19) 



T 



We can thus get 



aT = -^1^(0)^^(0)^0 = — - - 5/^7 , (6.20) 



as well as 



J-* 2,' 



Notice the relation 



^ = ^^JA I<;£)=b(")=0 = ^Bj^\g(0)^^m^^ , (6.22) 

which constitutes a valuable test in the numerical calculations. 

Finally, we find that the (non canonical) thermal conductivity is given by 

kT = —[e + p — 2fip — 26fi6p] + (Tap'^ + cfbS^^ + 2'yfi6fi , (6.23) 

CO 

where the term in the pressure p = e/2 has been added by hand to account for contact terms 
not directly implemented by the previous computations (see Herzog's review in [ll.J). 

The relations we have found between the different conductivities emerge quite naturally 
from the holographic setup. In the dual field theories they arise from Ward identity con- 
straints on the correlators (see, again, Herzog's review). 

Solving numerically the equations (6.5), (6.6) for the fluctuations of and and using 
the formulas above, we can calculate all the conductivities {aA,crB, i^,(yT, (3T) in terms 
of values of the dual gravity fields, 
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^^An alternative method consists in considering the hnear relations in (6.1) for dilFerent arbitrary choices 
of the values of the fluctuations at the horizon, in order to obtain enough equations to determine the various 
conductivities. Solving the system numerically we have checked that this method is stable w.r.t. those 
choices and it gives the same results as the method described above. 
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Figure 4: The real part of the electric conductivity for Sfi/fi = 0,0.8 (left plot, right plot) 
at Tc (dashed curves) and T > (solid curves). 



6.3 Conductivities in the normal phase 



As we have noticed in Section |4| following [H] , the normal phase of our model can be seen 
as a simplified holographic realization of a "forced" ferromagnet. Studying the conductivity 
matrix in this case is thus quite interesting, since it allows us to make some parallel with 
known results in ferromagnetic spintronics. 

In the left plot of Figure [i] we present the results for the real part of o"^ (the optical 
electric conductivity) as a function of the frequency of the external field perturbation for 
vanishing imbalance, 6fi = 0, i.e. the case considered in [15]. The different curves correspond 
to different temperatures; the dashed curve is at the critical temperature for the onset of 
superconductivity. At very large temperature, the conductivity is basically constant (quan- 
tities are normalized such that the constant is precisely equal to 1); this feature is peculiar 
of the gravity description we are using, 
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As the temperature is decreased, the conductivity is depleted at small frequencies. In 
fact, the imaginary part of cr^ (not shown) has a pole at w = 0; this is mapped by a 
Kramers- Kroning relation to a delta function for Re[(TA] at the same point (the solid line 
at u/T = in the figure). The delta function at zero frequency is due to the system 
translation invariance which, in charged media, causes an overall uniform acceleration and 
so an infinite DC conductivity!^ Since the area under the curves must be constant at different 
temperatures due to a Ferrell-Glover-Tinkham sum rule, the development of a delta function 
at cd = is compensated by a depletion of the conductivity at small frequencies 

In the right plot of Figure |4] we present the results for the real part of a a for no n- vanishing 



^^It depends on electro-magnetic duality of the four-dimensional Einstein-Maxwell theory on AdS4 |33| . 

^"This infinity is of course expected to transform into the standard Drude peak upon inclusion of impurities 
breaking translation invariance. Note that we are working with the fully backreactcd solution, so there is no 
dissipation as it happens in the probe limit, and as a result translational invariance is preserved. 
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Figure 5: The real part of the "spin-spin conductivity" ctb (left) and of the "spin conductiv- 
ity" 7 (right) for 6fi/fi = 0.8 at Tc (dashed curves) and T > Tc (sohd curves). 



imbalance, = 0.8. The general trend is clearly the same as in the balanced case. 

Nevertheless, the depth of the depletion at small frequencies is reduced with respect to the 
= case (remember that Tc is reduced too): the magnitude of the delta function is 
reduced, i.e. the DC conductivity is decreased by 

Concerning the real part of as (the "optical spin-spin conductivity"), it is exactly constant 
in absence of imbalance, 6fi = 0: there is no "net spin" in the system, so the conductivity 
is featureless]^ The unbalanced case, dj^/jj, = 0.8, is reported in the left plot of Figure [s] 
The similarity with the plot of a a, including the infinite DC conductivity, is transparent. In 



fact, it is clear (for example from equations (6.5), (6.6)) that in the normal phase (that is, 
at zero tp) the system enjoys the symmetry 

^^(/i, Sfi, u/T) = asiSfi, /i, u/T) , (6.24) 

which we also verified numerically with O{10^^) accuracy (at least). In particular, 

o'A = o'B for fi = 6^. (6.25) 

This is the same as stating that for perfectly polarized electrons (say all spins "up"), the 
electric conductivity cTcc (cta) equals the "spin-spin" one ags (o"_b)- In this case, then, we 
recover the zero-momentum result in ITOl. 



Actually, formula (6.24) is just a part of a set of more general relations among the vari- 



ous conductivities in the normal phase. In fact, it turns out that all the conductivities can 
be given once a single frequency dependent function / = f{u/T,6fi/fi) (a "mobility func- 



tion" for the charge/spin carriers) is known. Considering formulas (6.20), (6.21), (6.23), the 
^^The gravity vector Ba fluctuates on a black hole background charged under the other U{1)a- 
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conductivity matrix a in the normal phase turns out to b"'^ 

aT 



a 




(6.26) 



fP 



Pifp'^ + '^) -SpfpSp 
i^^Sp{f6p' + l)-pfp6p 



fp Sp 

ilE^Sp{f6p' + l)-pfp6p 
f6p' + 1 



These relations concern both the real and the imaginary parts of the conductivities. The 



matrix (6.26) reproduces in many instances the expectations in [8], for example despite 
their complicated explicit expressions, we have /3T = {5p/p)aT (since in the normal phase 
p5p = p6p, see Section |4]). The diagonal conductivities are related quadratically to the 
corresponding charge densities because the more the carriers are coupled the more they feel 
the external field, and, in addition, the bigger the charge density the stronger the transport. 
Again, the form (6.26) has been verified numerically with at least 0(10^'^) accuracy. 



The right plot of Figure [5] reports the results for 7, the mixed conductivity (the "spin 
conductivity" asc = o"cs), at 6p/p = 0.8. At 6p = this conductivity is identically zero: in 
absence of "net spin" , an external electric field does not cause the transport of "spin" (and 
the other way around). Instead, as Sp/p 7^ 0, we have the typical phenomenon at the roots 
of spintronics: an external electric field causes the transport of "spin" (and the other way 
around) with a non-trivial conductivity 7. 

Sp case (which from formulas (4.3), (4.4) is equivalent to p = 6p), 
1 



Note that in the p 



(6.26) implies 7 = cta — 1 = o"_b — 1: for perfectly polarized electrons the "spin conductivity" 
7 is equal to the electric conductivity a a, except for the constant contribution at large u/T 
(the "1" in formula (6.26) with our normalizations of the constant conductivity). The latter 
contribution marks the difference with the quasi-free electron case recently studied in [TO], 
where 7 = as at = Sp. The discrepancy is simply due to the fact that, differently from 
the quasi- free electron case, in the system at hand we have conduction even in the absence 
of net charge/spin density. 

In Figure |6] we present the real part of cr^ (left) and as (right) at fixed temperature for 
different values of 6p/p = 0, 0.8, 1.6 (dotted, dashed and sohd lines respectively),^^ Note the 



opposite behavior, dictated by formula (6.24), of the conductivities ctajCTb with increasing 



which determines the increase of the as DC conductivity with 6p/p. This behavior. 



^^We have left kT implicit just to save space: its expression, in terms of the other conductivities and of 
the thermodynaniical parameters in the normal phase (as given in Section [4]) can be immediately deduced 

de: 



from eq 

^■^We consider a generic setup where (5^//i is not constrained to be smaller or equal to one 



In QCD 

with up-down quark condensates, for example, Sfj, is given by the isospin chemical potential, while /i is the 
baryonic one, the two being in principle independent. 
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Figure 6: The real part of the electric conductivity a a (left) and of the "spin-spin conduc- 
tivity" as (right) for 5yu/yU = 0,0.8, 1.6 (dotted, dashed and solid lines respectively) at fixed 
temperature. 




(jjIT Cl>/T 

Figure 7: The real part of the "spin conductivity" 7 (left) and of the thermal conductiv- 
ity K (right) for Sfi/^ = 0,0.8,1.6 (dotted, dashed and sohd lines respectively) at fixed 
temperature. 



which has obvious physical origin in our system, is present in other contexts as well (see 
e.g. Pl]), where it is usually interpreted as a separation of the dynamics of charge and spin 
degrees of freedom. 

In Figure [7] we present the real part of 7 (left) and of the thermal conductivity k (right) 
at fixed temperature for different values of 6fi/ fi = 0, 0.8, 1.6 (dotted, dashed and solid lines 
respectively). Notice the non-monotonic behavior with 6fi/fi of both Re[7] and Re[K] in the 
small frequency region. Notice moreover that there is a delta function (whose coefficient is 
enhanced by 6fi) in the DC thermal conductivity due to momentum conservation (translation 
invariance) which forbids the relaxation of the heat current p!5]. This is refiected in a pole 
in Im[fi;] at a; = 0. 

Finally, in Figure [s] we report the real part of the thermo-electric conductivity aT (left) 
and of the "spin-electric" conductivity f3T (right) at fixed temperature for different values 
of Sfi/fi = 0,0.8, 1.6 (dotted, dashed and solid lines respectively). 
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Figure 8: The real part of the thermo-electric conductivity aT (left) and of the "spin-electric" 
conductivity (3T (right) for Sfi/fi = 0,0.8,1.6 (dotted, dashed and solid lines respectively) 
at fixed temperature. 
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Figure 9: The real part of the electric conductivity for = 0,0.8 (left plot, right plot) 

at Tc (dashed curves) and T < (solid curves). 



6.4 Conductivities in the superconducting phase 



The thermal behavior of the conductivities in the superconducting phase is similar to the 
one in the normal phase and it is shown in Figures |9l 



10 



In Figures [TT| [13} we present the results for cr^,(jB,7 in the superconducting phase for 
6fi/ fi = 0, 0.8, 1.6 (dotted, dashed and solid lines respectively) at constant temperature below 
Tc. The dotted line in Figure[ll|corresponds to the balanced case (5/i = 0) of [I5]. The optical 
electric conductivity at small temperature presents a pseudo-gap at small frequencies, while 
it relaxes to the normal phase value at large u. The corresponding imaginary part of the 



^^Since the system is really a superfluid, the spectrum is ungapped, containing the Goldstone boson of the 
breaking of U{1)a- This forbids the presence of a hard gap at T = [26]. Nevertheless, the conductivity 
is extremely small (almost exponentially) at small frequencies for sufficiently low temperatures, hence we 
speak about a "pseudo-gap" . 



28 




100 200 
co/T 




100 200 
co/T 



Figure 10: The real part of the "spin-spin conductivity" (Tb (left) and of the "spin conduc- 
tivity" 7 (right) for 6fi/fi = 0.8 at Tc (dashed curves) and T < Tc (solid curves). 
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Figure 11: The real part of the electric conductivity a a (left) and of the "spin-spin conduc- 
tivity" (Tb (right) for 6fi/^ = 0,0.8, 1.6 (dotted, dashed and solid lines respectively) at fixed 
temperature below Tc. 



conductivity (not shown) has a pole at w = 0, which translates in a delta function at the 
same point in Re[(TA] (the solid line at w/T = in the figure) due to a Kramers-Kroning 
relation. A part of this infinite DC conductivity is due to translation invariance and it is 



present also in the normal phase, as described in Section 6.3 Nevertheless, in Figure 12 (left 
plot) it is shown that the imaginary part of the conductivity has discontinuous derivative at 
Tc, and so the coefficient of the delta function has a jump across the phase transition. This 
means that a part of the delta function corresponds to genuine superconductivity and is due 
to the spontaneous breaking of U{1)a via the condensation of the charged operator dual to 

As is turned on and increased, the pseudo-gap in the real part of the electric con- 



ductivity Re[(jA] is reduced and eventually lost (dashed and sohd lines in Figure 11). This 
behavior has to be expected, since for large 6fi/fi the system at any fixed temperature 
will pass, with a second order, continuous transition, to the normal phase; accordingly, the 
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Figure 12: The discontinuity at in the imaginary part of the electric conductivity a a (left) 
and of the "spin-spin conductivity" cr^ (right), signaling a discontinuity in the magnitude of 
the delta function at w = in the respective DC conductivities. 



conductivity has a very similar behavior. Precisely the same pattern is seen in the opti- 
cal conductivities of some iron-based superconductors with increasing doping fraction (e.g. 
I35F 



Inversely, from the right plot in Figure 11 we see that the real part of the "spin-spin 



conductivity" Refcre] is more and more depleted at small frequencies as is increased. 



Also in this case it can be shown (right plot of Figure 12), that the imaginary part of 



the conductivity has discontinuous derivative at T,, the coefficient of the delta function at 



a; = presents therefore a jump across the phase transition temperature,^ That is, even 
if the condensing operator is uncharged with respect to U{1)bi due to the mixing with the 
U{1)a current in the superconducting phase, the as DC conductivity is enhanced. This 
strong coupling behavior is in stark contrast with the usual weak coupling picture, where 
the conductivity cxs is reduced by the decrease of the quasi-particle population due to the 
formation of the gap |36] . 



Figure 13 (left) shows that the "spin conductivity" 7 behaves qualitatively as as- For 
= it is constant and vanishing: in absence of "net spin" , an external electric field does 
not cause the transport of "spin" (and the other way around). Instead, a depletion at small 
frequencies opens as is increased and the DC conductivity becomes infinite. 

The "current polarization" P =Re[J^]/Re[J"^], measuring the difference of the two DC 
spin-polarized currents in absence of the "spin motive force" , is exactly equal to Re [7] /Re [a^] . 
In particular, it does not depend on the applied external field E"^. In the normal phase, the 



^^While in Fe-based superconductors the doping dependence does not directly translate into a chemical 
composition dependence as in cuprates, we chose to mention such an example due to the fact that Fe 
superconductors are believed to be s-wave (more precisely, S-| ), as the ones considered in this paper. 

^^The existence of this jump can be also verified, following section 4.2 of [TS], by considering the analytic 
form of the small frequency regime of the imaginary part of the conductivity slightly above the critical 
temperature, and checking numerically that the coefficient of the pole has a discontinuity when going slightly 
below Tc- The size of this discontinuity is of the same order of the one for a a- 
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Figure 13: The real part of the "spin conductivity" 7 (left) and of the thermal conductivity k 
(right) for 5yu/ fi = 0, 0.8, 1.6 (dotted, dashed and solid hues respectively) at fixed temperature 
below Tc. 



relations summarized in (6.26) imply that P = 6p/p exactly. In the superconducting phase. 



it can be checkeq^ that the qualitative behavior of P{Sp/p) is still increasing, but in a 
non-linear way (the power dependence being larger than one). For examples of current 
polarization behaviors in BCS systems, see e.g. [37]. The "optical current polarization", 
instead, can be shown to be non-monotonic with Sp/p. 

In the presence of the condensate, the symmetries encountered in the normal phase and 



encoded in formulas (6.24), (6.26) are explicitly broken (relations (6.20), (6.21) and (6.23) 



are preserved). In particular, (6.26) should be modified to account for the vev of the charged 
operator, which heavily infiuences the small frequency behavior of the conductivities. Un- 
fortunately, we did not find any simple modification that accounts for the numerical results. 

At this point, a comment is in order. It is clear that for non-extremal values of 6p/p, 
exemplified by the dashed lines in the figures of this section, the conductivities have a non- 
trivial (even non-monotonic) behavior for small frequencies]^ This is reasonable. In fact, the 
presence of the new energy scale provided by the imbalance can modify, with respect to the 
balanced case, the behavior of the system at frequencies related to 6p. We know from |33j 
that the behavior of the balanced system at zero charge density has just one regime at zero 
momentum: the hydrodynamic and the high frequency regimes coincide and the conductivity 
is constant. On the other hand, the condensate in the charged system changes drastically the 
behavior of the electric conductivity in the small frequency regime, which is now dominated 
by the superfiuid Goldstone mode and the pseudo-gap [HI |2S] . It is thus natural to expect 
a non-trivial infiuence of the extra IR scale 6p on the transport properties of the theory. In 



^^Since the DC conductivities are infinite in the superconducting phase, we estimate the ratio Re[7]/Re[o'^] 
by evaluating the corresponding ratio of imaginary parts at very small uj /T, relying on the Kramers-Kronig 
relation for the result to translate in the respective behavior of the real parts. 

^^We think that this behavior is not due to numerical effects because we checked it with two independent 
codes and two different methods for the calculations. 
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Figure 14: The real part of the thermo-electric and "thermo-spin" conductivities (left plot, 
right plot) for = 0,0.8, 1.6 (dotted, dashed and solid lines respectively) at fixed tem- 

perature below Tc. 



particular, from the plots of this section it seems that the DC electric conductivity is further 
enhanced by 5/i w.r.t. the extrapolated behavior of the normal phase, possibly due to a 
second light mode. 

The thermo-electric, "spin-electric" and thermal conductivity follow from the plots above 
and formulas (6.20), (6.21), (6.23). The corresponding plots are reported in Figure [l3| (right) 
andlH 

Finally, in Figure 15 we report the plot of the "pseudo-gap" frequency Wgap below which 
the real part of the optical electric conductivity a a is essentially vanishing,^ at constant tem- 
perature but increasing The behavior of uj^s_^{5iJi/ jj) is clearly non-linearly decreasing. 
This has to be compared to the case of ordinary unbalanced superconductors, where the gap 
A is constant in 5/i at T = 0. In [15] it was pointed out that Co'gap(5/i = ^)/T^ ~ 8 as in 
some measures in high superconductors. We find that a;gap(5/i//i)/Tc is not approximately 
constant, but a decreasing function of which substantially deviates from the value 8 
even in the range of where it is still reasonably well defined. 



7 Comments on possible string embeddings 

If the results we have found on "charge" and "spin" transport properties have some degree 
of universality, as we have argued above, the same is not necessarily true for equilibrium 
properties. The (T, 5^) phase diagram, for example, can vary for different choices of the 
scalar potential, as it happens in the balanced case (see for example [38]). This variation 
accounts for different microscopic properties of the dual field theories which are however 
unknown. The latter could be unvealed only by embedding the bottom-up models in full- 

^^To be precise, we used the numerical threshold value Re[crA(wgap)] = 0.005. 
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Figure 15: The "pseudo-gap" frequency function of at fixed T. The error 

bars are associated to numerical uncertainty and the plotted line emerges from a fit with a 
quartic polinomial. 

fledged string or M-theory constructions. 

The embeddings depend both on the spacetime dimensionality and on the microscopical 
details of the dual fleld theory - essentially on the kind of U{1) gauge flelds entering the game 
and on the precise nature of the order parameter. The latter can in fact be a condensate of, 
say, adjoint (or more generic two-index representations of some gauge group) or fundamental 
fermionic degrees of freedom and the string or M-theory embedding would strongly depend 
on whether one or another possibility is realized. 

7.1 Adjoint fermion condensates 

In supersymmetric contexts, this is a case where, say, gluino bilinears break some f/(l)i? 
symmetry of the theory]^ In this case we should try to see whether our minimal 3 -|- 
1 dimensional gravity model can arise from a consistent Kaluza- Klein truncation of lld- 
supergravity on some compact seven-manifold, in the same way as it happens (see [101 SI] ) for 
the balanced model introduced in [H]. Isometries of the seven-manifold, in fact, are mapped 
into global (R-) symmetries of the dual field theory. At least within known consistent KK 
truncations (see e.g. J41j) it seems that an embedding of a fully backreacted model containing 
the same fields as ours is possible only provided at least another non-trivial real scalar field 
is present. This is true also in the normal phase, i.e. in the absence of our complex scalar 
field. 

■^°J7(l)fl superfluidity driven by gluino Cooper pairs in A/" = 4 SYM at finite density has been first 
considered in [55] , 
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To see this explicitly, let us consider the KK reduction of lld-supergravity on a seven 
sphere 5*^. This gives rise to gauged d = 4, A/" = 8 5*0(8) supergravity, which can be further 
truncated to a gauged Af = 2 model where the bosonic sector consists of the metric, four 
Maxwell fields A^, three '"dilatons" (pi and three "axions" [l^. The resulting Lagrangian 
density reads 

{y^g)-'C = R-li2 + %'(cosh0.)] - J E E e'^'-'^^iPy + . . . , (7.1) 

i=l A=l i=l 

where g is the coupling, a\ are constants and we have not included the contribution of the 
axions. This action cannot be further truncated to a model with constant dilatons and two 
independent non-trivial Maxwell fields. Such kind of truncation is only allowed if the U{1) 
fields are identified up to some constant (hence in the balanced case), and this constraint 
applies also to more general KK truncations 

The same conclusion holds in the d = 5 case, i.e. considering abelian ^2] or non-abelian 
gauged supergravities |13] obtained from consistent KK truncations of JIB supergravity. For 
example in the first reference in [H] it is shown how a f/(l) ^-charged RN-AfiS" black hole 
solution like the one describing the normal phase in our model, can be embedded within 
Romans' J\f = 4 SU{2) x U{1) gauged supergravity [13] (and thus in type IIB), only if the 
two Maxwell fields (and so the corresponding chemical potentials) are identified modulo a 
constant. Again, a setup where the two fields are independent can only be realized adding 
a non-trivial real scalar (dilaton-like) field. Embeddings of balanced superconductors |14] , 
instead, do not require this extra scalar to be present. 

7.2 Fundamental fermion condensates 

^ Matter fields transforming in the fundamental representation are introduced in the holo- 
graphic correspondence by means of flavor D-branes. Models of holographic p-wave super- 
conductivity have been actually embedded in probe flavor brane setups US] • Here we 
want to focus on the s-wave case. Having in mind 4d QCD-like models, one could con- 
sider, say, a non-critical 5d string model with Nc D3 and Nf spacetime filling D4-anti-D4 
branes [IHl I1Z|. The low-energy modes of the D3-branes would be the SU{Nc) gluons, and 
the D4-anti-D4 branes would provide the left and right handed fundamental flavor fields 
(the quarks). The model contains a complex scalar field (the would-be tachyon of the open 
string stretching between branes and anti-branes) transforming in the (anti)fundamental of 
SU{Nf)L X SU{Nf)R: its condensation (eventually triggered by a running dilaton which can 
account for confinement) drives the breaking of the chiral symmetry down to SU{Nf) and 
thus it is dual to the chiral condensate of fundamental fermions |1H1 SSI HZ] . The model (at 

^"'^We are grateful to Erailiano Imeroni for his contributions to this part of the project. 
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least a simplified version of it can provide AdS^ flavored solutions (at zero temperature 
and densities) with trivial tachyon and constant dilaton. These are possibly dual to phases 
of the theory in Banks- Zacks-like conformal windows. The corresponding finite temperature 
versions have been studied in [i^ . 

Let us start considering the Nf = 1 case. In addition to the above mentioned complex 
scalar field r, the DBI brane-antibrane action [^01 EI] contains two Maxwell fields and it 
is coupled with the dilaton. The scalar field r is charged under a combination of the two 
gauge fields (the "chiral" U{1)a) and uncharged under the orthogonal combination (the 
"baryonic" U{1)b)- The chiral symmetry is actually anomalous and this is accounted for 
by other terms in the D-brane action [17]. One can consider Nf > 1 setups, e.g. Nf = 2 
ones, enhancing the preserved flavor symmetry to an SU{2) x U{1)b one. Starting from 
the above mentioned flavored- AciS" solutions one could turn on a chemical potential for a 
U(l)i C SU{Nf) "isospin" fleld as well as a baryonic one. Analogously to our condensed 
matter setup, one could then study the stability of the system under fluctuations of the 
complex scalar fleld r. The latter would be dual to, say, a ud mesonic condensate and thus 
the role of 6fi and /i would be played by /i^ and /i/ respectively. Let us focus on the Nf = 1 
case for simplicity, treating the axial U{1)a as if it were be a genuine symmetry replacing 
U{l)j, i.e. neglecting all the terms in the brane action which are related to its anomaly. 

The Dp-brane-anti-brane DBI action in string frame reads [501 151] 

S = -T„ 



j dP+^x e-*V^(|r|) [V-detA(^) + V-detA(«)] , (7.2) 



where Tp is the D-brane tension, $ is the dilaton and ^(|t|) is the open string tachyon 
potential. Let us assume that a closed string tachyon, which could be present in non-critical 



Type setups, is eventually frozen out. The matrices A*^^)'*^^) in (7.2) are deflned as 

^MN = P[9 + B2]mn + 2WF«jv + na'iDMrriD^r) + Txa\D^Ty{DMr) , 

= OmA^S - QnA^S , Dmt = {dM + %A^S - V , (7.3) 

where i = L,R label the brane or antibrane, P[«] denotes the pullback on the D-brane 
worldvolume, g is the metric and B2 is the NSNS antisymmetric two-form. Notice that, as 
anticipated, the scalar fleld is charged only under the axial U{1)a_ combination 

Am = - ^1? • (7.4) 
It is instead uncharged under the barionic U{1)b combination 

Bm = A^S + . (7.5) 

In the case of brane-antibrane pairs in flat spacetime, string fleld theory gives a tachyon 
potential of the form 

y(|r|) = e™'"^'l"l' , withm2 = -^. (7.6) 



35 



This expression could be affected by non-trivial field redefinitions and by the fact that the 
branes have to be put on curved spacetimes (see e.g. a discussion in [47J). Anyway, we will 
take this expression as a guideline. 

As for the embedding described in the previous subsection, when independent Maxwell 
fields are both turned on, it seems not possible to have solutions with non-trivial dilaton (it 
is understood that the whole gravity action, say for p = 4, will be a "bulk+brane" one with 
the standard kinetic term for the dilaton). Notice that in the model the Maxwell fields are 
coupled by the non-linear structure of the DBI action and thus spintronics effects could be 
present also in the non-fully backreacted case. 

Let us just notice, finally, that on a closed string background where B2 = and the metric 
is diagonal, the low energy effective Lagrangian density coming from (7.2) at the quadratic 
level in the fields is (see eq. (10) in the second paper in [51] and use our redefinitions (7.4), 
(Q) 

C ^ -Tp(2W)e-*v^ [-Tra'F^ - na'Y^ + \Dr\^ + mVl'] , (7-7) 

where F = dA, Y = dB and Dt = {d — iA)T. This has, modulo coefficients and the overall 
dilaton coupling, the same form as the matter part of our gravity model. 
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A The Chandrasekhar-Clogston bound 

Let us consider a Fermi mixture with two species u and d (e.g. dressed electrons of spin up 
and down in metallic superconductors) with different chemical potentials fiu and /x^. Let us 
define the mean chemical potential /x and the chemical potential imbalance as 

= \{l^u + fJ'd) Sfi = }-{nu- fJ^d) ■ (A. 8) 
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Assuming analiticity, in the grand-canonical ensemble, the Gibbs free energy at zero 

temperature can be Taylor expanded for 5/i ^ yU as 

Q{6fi) = Q{0) + Q{0y6fi + ^Q"{0)6fi^ + 0{6fi^) . (A.9) 

The first (resp. second) derivative of Q w.r.t. S/j, defines the population imbalance 6n (resp. 
the susceptibihty imbalance '^^ 6x) 

Sn^n.-n, = -—, = ^ = • (A-IO) 

In BCS theory, the normal phase at T = and 5^ <^ ^ has a population imbalance given 
by 5njq ~ ppSfJ^, where pp is the two-Fermion density of states at the mean Fermi surface 
E = Ep = p- From this expression we easily find that, in the normal phase, the free energy 
expansion (A.9) reduces at leading order to 

n^{6p)^nNio)-^PF6p\ (A.ii) 

At T = 0, the homogeneous BCS superconducting phase, characterized by Cooper pairs of 
zero total momentum, has an equal number of particles of species 1 and 2: Sns = 0. Thus 
the free energy just expands as 

ns{Sp)^ns{o). (A.12) 

The difference between the two free energies is given by 

n^i^p) - ns{6p) ^ n^io) - ns{o) - \pf^p'. (a.i3) 

Now, using the standard BCS result, fiAr(O) — ^^^(O) = pi?AQ/4, where Aq is the gap param- 
eter at T = 0, it follows that 

n^i^p) - ns{6p) ^ ^-PF^I - ^PfSp^. (A.14) 

This shows that at T = the superconducting phase is favored (i.e. its free energy is less 
than the free energy of the normal phase) only if the Chandrasekhar-Clogston bound 

Sp < Spi, Spi = , (A. 15) 

v2 

is satisfied. 



■^^In the case where the chemical potential imbalance is induced by the Zeeman coupling with a magnetic 
field, this is actually the magnetic susceptibility. 
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B Equations of motion in d + 1 bulk spacetime dimen- 
sions 

Let us consider the generalization of our model to d+ 1-dimensions. The action reads 
1 



9k^ 



9 



n + ^^^^^ - \f,,f^' - ^Y,,Y^^ - vm - lav^ - ig^vr 



(B.l) 



The ansatz for the spacetime metric is 



ds" = -g{r)e-^^'-^de + —dx^ + — - , (B.2) 

g{r) 

together with an homogeneous ansatz for the fields 

V' = i^{r) , Aadx"" = (f){r)dt , Badx" = v{r)dt . (B.3) 
The equation of motion for the scalar field reads 

\g r 2J 2 g g^ 
Maxwell's equation for the field gives 

0"+0'fi^ + ^)-2^^ = O. (B.5) 
\ r 2 / g 

Einstein's equations reduce to 

1 e>^{4^'^ + v'^) (d - 1) g' l (d-l)(d-2) _ d{d - 1) m , q'r<t>^e>c ^ . 
2^ Ag 2 gr 2 2gL'^ 2g 2g^ 

Finally, Maxwell's equations for the additional gauge field read 

y" + y'f{^^ + ^)^0_ (B.7) 
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B.l The normal phase 

The gravity solution corresponding to the normal phase in the dual d-dimensional field theory 
is the U{iy - charged Reissner-Nordstrom (RN)-A(i5'(i_|_i black hole 



9{r) 
At 
Bt 



-g{r)dt' + —dx^ + 
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(B.8) 
(B.9) 
(B.IO) 
(B.ll) 



where r = r^^ is the position of the outer horizon. 

The charge densities of the dual field theory are related to the subleading behavior of the 
bulk Maxwell fields as 



P 



1 {d-2)i,r % 
The black hole temperature is 



d-2 



5p 



1 {d - 2)5/ir^" 



d-2 



ld-1 



(B.12) 



d 



{d-2y {fi^ + 6fi^)L^ 



2r2 



(B.13) 



The Gibbs free energy density (hence the pressure) is given by 



1 



cu = —p 



H 



id-2) {fi^ + 6fi^)L^ 
2{d-l) 



H 



(B.14) 



Consistently, the charge densities in (B.12) are obtained as 

du 



P 



du 



Sp 



dSp 



(B.15) 



The energy density is given by 



d-l rjj 



{d-2) (/i^ + V)^' 



1 + 



H 



(B.16) 



This satisfies the relation e = [d — l)p related to the vanishing of the trace of the stress 
energy tensor. 
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B.1.1 Near horizon geometry 



The Reissner-Nordstrom geometry is interesting as T — t- 0. In this limit the horizon radius 
has a fixed value at 

To find the near horizon metric take the series Taylor expansion of the blackening factor 

g{r) ~ g{rH) + g'{rH)f + ^g"{rH)f'^ , (B.18) 
where again r = rn + r with f — > 0. We find that 

9{rH) = 0, g'{rH)r^T = 0, /(r^) = M^zll . (b.19) 
The near horizon metric reads then 

dsl^r horizon " "^(^ " ^)^dt^ + ^rff^ + ^(^_^)-2 ^^' ' (^-20) 

from which we recognize the AdS2 x i?*^"^ metric. The AdS2 radius squared is L^g) = 
LV(rf(rf-l)). 

B.2 Criterion for instability 

Let us consider the stability of the above solution at T = under fiuctuations of the charged 



scalar field. The equation one has to consider is given in (B.4). The background is the 
extremal doubly charged KN-AdSd+i- In the near horizon limit it is easy to show that the 
equation of motion for reduces to that of a scalar field of effective mass 



"^(2) = - ' x = 5fi/fi, {B.21] 

\ ) 1 + X 



on an AdS2 background of radius L^2) = L'^/{d{d — 1)). 

The background is unstable in the limit if the BF bound in AdS2 is violated, i.e. if 
L^2)"^(2) < —1/4. This is equivalent to the condition 



which generalizes our formula (4.14) valid for d = 3 
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